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ABSTRACT 



We provide hot star wind models with radiative force calculated using the solution of comoving frame (CMF) radiative 
transfer equation. The wind models are calculated for the first stars, O stars, and the central stars of planetary nebulae. 
We show that without line overlaps and with solely thermal line broadening the pure Sobolev approximation provides 
a reliable estimate of the radiative force even close to the wind sonic point. Consequently, models with the Sobolev line 
force provide good approximations to solutions obtained with non-Sobolev transfer. Taking line overlaps into account, 
the radiative force becomes slightly lower, leading to a decrease in the wind mass- loss rate by roughly 40%. Below the 
sonic point, the CMF line force is significantly lower than the Sobolev one. In the case of pure thermal broadening, this 
does not influence the mass- loss rate, as the wind mass-loss rate is set in the supersonic part of the wind. However, when 
additional line broadening is present (e.g., the turbulent one) the region of low CMF line force may extend outwards 
to the regions where the mass-loss rate is set. This results in a decrease in the wind mass-loss rate. This effect can at 
least partly explain the low wind mass-loss rates derived from some observational analyses of luminous O stars. 

Key words, stars: winds, outflows - stars: mass-loss - stars: early-type - hydrodynamics - radiative transfer 
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One of the most important galactic populations consists 
of massive stars, because these stars dominate the spec- 
tra of many galaxies and contribute significantly to the 
mass and momentum input into the interstellar matter. 
Moreover, massive stars end their active lives in gigantic 
explosions such as supernovae or even possibly as the pro- 
gcnitors of ganima-ra y bursts (see Wooslev & Hcger 200g, 
lYoon fc Langeij [20051 ). producing huge amounts of heavier 
elements. 

An important property of hot stars that significantly 
influences the i r final stages is the stel l ar wind (see, e.g. , 
lOwockil [200I iKrticka fc KubatI [2007aL iPuls et all [2008bl 
for reviews dedicated to hot star winds). However, in stel- 
lar evolution calculations it is usually unnecessary to know 
detailed wind properties, but just the amount of mass ex- 
pelled from the star per unit of time (mass-loss rate) as a 
function of stellar parameters (e.g., mass, effective temper- 
ature, radius, surface metallicity) . However, for many hot 
stars we simply cannot estimate their true mass-loss rate 
with the precision necessary to calculate evolutionary mod- 
els. The situation may be less problematic for luminous O 
stars, for which relatively good agreement between theo- 
retical predictions and observational results seems to exist 
(iPauldrach et al.ll2001l JVink et al.|[200l IKrticka fc Kubatj 
|2004 hereafter iPaper t ). 

However, the agreement between theoretically predicted 
mass-loss rates and those derived from observations may 
be an illusion caused by the neglect of some physical ef- 
fects in the wind, such as clumping ([Bouret et al.l 120031 
iMartins et al.l 120051) . As a result, the true mass- loss rates 



of O stars may be a few times lower than the standard wind 
theory pr edicts. This seems to b e supported by the obser- 
vations of iFullerton et all (120061 ) of weak wind line profiles 
of P V. Last but not least, the unexpected occurrence of 
symmetrical X-ray line p rofiles seems to requ ire relatively 
low wind mass-loss rates (jKramer et al.l[2003l) . 

Possibly unreliable estimates of hot star wind mass-loss 
rates are also problematic because altough more realistic 
evolutionary stellar models can be calculated, by includ- 
ing, e.g., rotation and magnetic fields, the wind mass- loss 
rates remain uncertain. Ideally, all observational indica- 
tors of mass-loss rate and theoretical models should find 
and predict similar mass-loss rates. From the observational 
point of view, more detailed models of line formation in 
inhomogeneous media may be necessary to obtain reliable 
line profiles, and conseque ntly also e stimate mass-loss rates 
(Oskinova et al. 2007, Sundqvist et al. 2010) . 

From the theoretical point of view, disagreement be- 
tween theory and observations would imply that some of 
the assumptions used for the hot-star wind modeling are 
inaccurate. Part of the disa greement may be cause d by 
using incorrect abundances ([Krticka fc KubatI l2007b[) . al- 
though the reason for a disagreement remains mainly un- 
clear. A thorough inspection of all assumptions involved 
in the modeling is therefore strongly needed. As a first 
step in this direction, we studied the influence of X-rays 
on the wind structure of hot stars. It seems that X-rays 
alone cannot entirely e xplain the disagreement between the- 
ory and observations (jKrticka fc Kubati 12009 ^ as their in- 
fluence on wind mass-loss rates is small and they do not 
strongly affect the ionization fraction of many important 
ions, especially that of P V. On the other hand, the modified 



J. Krticka and J. Kubat: CMF models of hot star winds 



ionization equilibrium may affect the X-ray l ine formation 
(jOs kinova e t al.ll2006Llkrticka fc Kubat "2009V and too ling 
cooling time in the pos t-shock region (Cohen ct al. 2008, 
iKrticka fc KubatI 2009D may cause the so-called "weak 
wind problem" (iBouret et all [200l iMartins et al] |2004 
iMarcolinoet al.ll2009D . 

One of the most important approximations in the hot- 
star wind modelin g is the Sobolev approximation (jSobolevI 
Il947l ICastoiJll974 ). which enables us to solve the line radi- 
ation transfer analytically. Some studies confirm its appli- 
cability in the supersonic part of smooth line-dr iven winds 
()Hamannlll98TllPauldrach et al.lll986llPulslll987t) . However, 
the applicability of the Sobolev approximation is question- 
able especially in the regions close to the photosphere ba- 
cause of the existence of strong source function gradients 
(jOwocki fc P uis 1999). On the other hand, some models 
avoid using the Sobolev approximation and use only the 
comoving-frame (hereafter CM F) method of solving the ra - 
diative transfer equation (e.g.. lGrafener fc Hamannll2005D . 

We decided to test the applicability of the Sobolev ap- 
proximation and include the CMF solution of the radiative 
transfer equation in our wind models. In this first paper of 
a series, we describe our method, and study the applicabil- 
ity of the Sobolev approximation using models neglecting 
continuum opacity. 



2. Basic model assumptions 

The models used in this paper are based on the NLTE wind 
models of Krticka fc Kubat (2004, hereafter iPaper J ). Here 
we summarise only their basic features and describe the 
inclusion of CMF line force. 

We assume a spherically symmetric stationary stellar 
wind. The excitation and ionization state of the consid- 
ered elements is derived from the statistical equilibrium 
(NLTE) equations. Ionic models are either adopted 
from the TLUSTY gr i d of model stellar atmospheres 
(|Lanz fc Hubenvl l2003l |2007[ ) or are prepared by us 
using the data from the Opac it y and Iron Projects 
(ISeatonI I1987L iFernlev et all I1987L (Luo fc PradhanI Il989l 



Sawev fc Berringtonlll992l ISIeaton et al.lll992llButler et alj 
19931 iNaha r fc Pradhan '1993'. 'Hummer ct al.' '1993'. 
Bautistal 119 96. Nalia r fc P ra dhan 1996, Zhan g 1996, 
Bautista fc Pra dha^ Il997l IZhang fc PradhanI ' Il997l 



19991) 



Chen fc PradhanI I1999D . As in IPaper J . the solution of 
the radiative transfer equation for NLTE equations is 
artificially split into two parts, namely the radiative 
transfer in either the continuum or in lines. The solution 
to the radiative transfer equation in continuum is based 
on the Feautrier method in the spher ical coordinates 
(jMihalas fc Hummeil 11974 iKubatI 1 19931) . and the line 
radia tive transfer i s solved in the Sobolev approxima- 
tion (jCastoiJ 11974 iRvbicki fc Hummeil |1978[) neglecting 
continuum opacity and line overlaps. 

In contrast to our previous models, the radiative trans- 
fer in lines used for the calculation of the radiative force is 
solved in the CMF (see Sect. |2|) neglecting the continuum 
opacity. The line radiative force is calculated directly from 
the true chemical composition, NLTE ionization and exci- 
tation bal ance, and CMF flux using data from the VALD 
database (jPiskunov et al.l 1995. Kup ka ct al. 19 93). We do 
not use the line-strength distribution function parameter- 
ized by force multipliers k, a, and S. 



The flux at the surface (used as the lower boundary 
condition for the radiative transfer in the wind) is taken 
from the H-He spherically sy mmetric NLTE model stellar 
atmospheres of lKubatI (|2003l and references therein). 



3. CMF calculation of the radiative force 

The radiative force is calculated using the solution of 
the spherically symmetric CMF radiative transfer equation 
(lMihalaslll978l . Eq. (14.99)) 

dl[v,^,r) 1 — ^^ 57(j/, /x, r) 
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Vr dr 
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= 'q{v,r)-x{v,r)I{v,ii,r), (1) 



where /(i^, /j,, r) is the intensity seen by the observer moving 
with the wind at the radial velocity Vr, i^ is the frequency, 
/i = cos 9, 9 is the direction between the given ray and 
the radial direction, and 77(1/, r) and x('^i '') ^'^^ the line 
emissivity and opacity, respectively, given by 
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where rii and Uj are the number densities of individual 
states with statistical weights gi and gj corresponding to 
the line transition i ^ j with oscillator strength /^ and the 
line-profile (pij{v), and me is the electron mass. Assuming 
the line profile to be Gaussian produced by thermal broad- 
ening only, fij{v) is given by 
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where Vij is the laboratory line frequency, and the line 
broadening is given by 



Aiy„ = ^ 



2fcT 



(4) 



and ma is the mass of a given atom. The number densities 
of individual levels in Eqs. ^ are calculated from statistical 
equilibrium (NLTE) equations. 

Writing Eq. ((H, we neglected advection and aberration 
terms, which is justifiable in non-relativistic flows (see, e.g., 
iKorcakova fc Kubltll2003l ). We also note that by neglecting 
the spatial derivatives o f intensity in Eq . ([Ij we obtain the 
Sobolev approximation (Castor 2004) . 

Following Mihalas et al. (1975 ) , we rewrite Eq. ([1]) for 
rays with an impact parameter p 



± 



dz 



cr 



dWr 



Vr dr 



dv 



l-/i^ 

= r]{v,r)-x{v,r)I^{v,p,z), (5) 



dl^{v,p,z) 



where -I- and — refers to radiation flowing toward and away 
from the observer, respectively, r = {p^ + z'^)^''^, and z 
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is the distance along the ray. We transform Eq. ([5]) using 
intensity-like and flux-like variables 

u{iy, p,z)^ - [/+ {v, p, z) + r {v, p, z)] , (6a) 

v{iy,p,z) = - [l+{iy,p,z)- r{i^,p,z)] , (6b) 

to obtain a system of partial differential equations 

1 du{iy,p,z) dv{v,p,z) 

- l{v,P, z) = -v{v,p, z), 



approach. We calculate the ratio of the CMF and Sobolev 
line forces 



X(i^, r) dz 
1 dv{h',p,z) 



diy 



x(i^, J') dz 
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The system of equations in Eq. (O is s olved numeri- 
cally u sing the long characteristic method of iMihalas et al.l 
(jl975|) . which we modified slightly for the present purpose 
(see Appendix R)) . As we are interested in the calculation 
of the radiative force us hig the v variable at a particular 
grid point, in contrast to lMihalas et al.l (|l975l ) we specify v 
at grid points, and u in the middle between them. 

In our numerical solution of Eq. ([7|), we use the 
same spatial grid as for the solution of hydrodynamical 
equations. The spacing of the frequency grid is Ai^d = 
V ^J {2kTQ / mc) / {cfi)), where Tq is the pre-specified ex- 
pected minimum wind temperature, mc is the atomic mass 
of artificial metallic atom, and /d is the multiplicative fac- 
tor (see below). The CMF radiative transfer equation is 
solved only for selected frequencies from the frequency grid 
that lie close to some line. The selection of frequencies is 
controlled by two in teger numbers ud, and NCERV (cf. 
iHillier fc Milleil[l998f) . For each line, we select frequencies 
that lie within tid line Doppler widths Avij. Redward of 
the center of each line, we select each NCERV frequency 
up to the frequency corresponding to the Doppler shift for 
the wind terminal velocity. The numerical test showed that 
a sufficiently precise value of the radiative force can be de- 
rived for the value of parameters mc — 60r7iH, where mn is 
the mass of hydrogen atom, /d = 2, no = 5, NCERV=30, 
and typically Tc = 10 000 - 20 000 K. 

The radiative force is calculated as an integral 
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As the calculation of the CMF radiative force is rather 
time-consuming, we do not calculate /™^ during each it- 
eration of hydrodynamical variables, but adopt a different 
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By the Sobolev line force ff°^^ we mean here the force 
calculated by assuming the Sobolev approximation for ra- 
diative transfer, neglecting line overlaps and using true 
line opacities and t he emer gent flux from the underlying 
stellar atmosphere (jPaoer J . Eq. (25) therein). Unless the 
base density is known with a precision better than about 
30%, we calculate c*-''^^ only when the estimate of the base 
density is changed, and keep c^^^ fixed during the subse- 
quent iterations of the hydrodynamical structure. When the 
base density is known with a higher precision, we calculate 
^CMF g^j^gj. each change of the hydrodynamical structure. 
Moreover, because we solve the hydrodynamical equations 
using the Newton-Raphson method, we have to calculate 
the derivatives of f^^^ with respect to individual hydro- 
dynamical variables. These derivatives are approximated 
using the derivatives of the Sobolev line force ff°^ multi- 
plie d by c^'^ ^. The force term in the critical point condition 
fsee iPaperB ) is also multiplied by c^^^ . 

We note that direct use of Eq. (fT3)) causes instability in 
the model convergence. The reason for these convergence 
problems may be numerical, but this b ehavior may also be 
connected with line-d riven instability (jOwocki et al.lll988L 
iFeldmeier et al.lll997[ ). To avoid this (since we are seeking 
stationary solution and not evolution with time) we intro- 
duced a weak smoothing of c*"''^^. 



^CMF 
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2c; 



CMF 



^CMF , CMF\ 



(14) 



where c^'^^ is the value of c*^^^ at a given grid point d (as 
for d — 1 and d + 1) and we use cj'^^ instead of c™^ in 
the models. Our numerical tests showed that the smoothing 
Eq. ([T4| does not significantly affect the resulting radiative 
force. 



4. Studied model stars 

In our study, we selected three types of stars to study more 
carefully the Sobolev approximation in different wind envi- 
ronments (see Table [ij. 

The stellar parameters of the first stars were ob- 
tained according to an evolutionary calculation of initially 
zer o-metallici t y sta r with initial mass 50 M© derived by 
M arigo et al.l ( 200 it) . For these models, we assumed a stel- 
lar wind driven purely by CNO elements (which appear on 
the stellar surface due to mixing) with a mass-fraction of 
CNO Z = 10-3. 

The stellar parameters (effective temperatures and 
radii) of an O star sample were der ived using the model 
atmospheres with li i ie blanketing (iRepolust et al.l |2004 
iMarkova et~aI1 12004 iMartins etall I2005D . Stellar masses 
were o btained using e volutionary tracks eith er by ourselves 
using ISchaller et al.l |1992| tracks) or by IMartins et al.l 
2005f). For these stars, we assumed a solar chemical com- 



position (jAsplund et al.ll2005D . 

The stellar paramete rs of central sta rs of planetary 
nebulae were taken from iPauldrach et al.l (J2004t ). who de- 
rived them fr om UV spectros c opy. H elium abundance was 
adopted from lKudritzki et al.l (|l997l ). for other elements we 
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Table 1. Radius i?, , mass M, and the effective temperature 
Toff of studied model stars 



Star (model) 


[Rq] 


M 

[Mq] 


Tcff 

[K] 


First stars 


M500-1 
M500-2 
M500-3 
M500-4 


11.1 
33.7 
72.0 
303 


50 
50 
50 
50 


50 000 
29 900 
20 600 
10100 


stars 



^ Per 

t Ori 

15 Mon 



COph 
68 Cyg 
19 Cep 



HD 24912 
HD 37043 
HD 47839 
HD 54662 
HD 93204 
HD 149757 
HD 203064 
HD 209975 



14.0 
21.6 
9.9 
11.9 
11.9 
8.9 
15.7 
22.9 



36 
41 
32 
38 
41 
21 
38 
47 



35 000 
31400 

37 500 

38 600 
40 000 
32 000 
34 500 
32 000 



Central stars of planetary nebulae 



NGC 2392 

NGC 3242 

IC 4637 

IC 4593 

He 2-108 

IC418 

Tc 1 

NGC 6826 



1.5 
0.3 

0.8 
2.2 
2.7 
2.7 
3.0 
2.2 



0.41 
0.53 
0.87 
1.11 
1.33 
1.33 
1.37 
1.40 



40 000 
75 000 
55 000 
40 000 
39 000 
39 000 
35 000 
44 000 




0.5 
0.4 



HD 24912 

HD 93204 

HD 93204, S=S„ 



0.1 



10 



100 



v,/, 



Fig. 1. The ratio of the CMF to Sobolev line forces given 
by Eq. (J13p as a function of the wind velocity plotted in the 
terms of the sound speed for two selected stars. The dashed- 
dotted line denotes a model with constant source function 
and constant level populations equal to their values at the 
sonic point. Arrows indicate the thermal speed of selected 
ions. 



assum ed a solar chemical composition (after lAsplund et al.l 
I20p5[) ■ which was for some stars slightly modified according 
to lPauldrach et al.l ()2004D . 



5. Comparison of CMF and Sobolev wind models 

We calculated wind models with both CMF and Sobolev 
line forces and compared the final wind structure. The re- 
sulting ratio of the CMF to Sobolev line forces c'~^^^ for 
selected stars is shown in Fig. [T] We note that the Sobolev 
force was calculated using the flux from the stellar atmo- 
sphere and by neglecting line overlaps. 



For very low wind velocities Vr < 0.1a (where a^ — 
2fcr/mH), the CMF force is large, c^^^ > 1. This is 
most likely partly connected with the boundary conditions, 
which are not completely c ompatible with the wind (cf., 
iNoerdlinger fc Rvbickilll974[ ). 

For velocities of about one tenth of the sound speed, 
there is an apparent minimum of c . In some cases, the 
ratio c could even be negative, which corresponds to 

a negative CMF radiative force. The Sobolev approxima- 
tion is not applicable in this region, but a low value of the 
radiative force is also connected with positive source func- 
tion gradients. For subsonic velocities, the Doppler shift is 
less important, and the line radiative transfer is given ba- 
sically by the static radiative transfer equation. In the op- 
tically thick regions, for frequencies corresponding to line 
transitions it follows from Eq. (j7bp u Ri S, from Eq. ([7a| 
V ~ — (l/x)dS'/dz, and the radiative force is proportional 
to the negative of the derivative of the source function 
fgy^ ~ -dS/dz (see Eq. dH]), and INoerdlinger fc RvbicO 
|1974[ ). Because the line source function increases here (see 
Fig. [2]), the line radiative force at low velocities may even 
be negative. For a constant source function, the mini- 
mum of c''-^^^ close to the star is significantly weaker (see 
Fig. [1]). The source function minimum below the sonic 
point is caused by a local temperature minimum, because 
the line source function of optically thick lines (which are, 
consequently, in detailed radiative balance) close to the 
star 5* « rij/ni ^ {rij/m)* (asterisk denotes LTE value) 
depends on temperature. Another source function mini- 
mum for non-Sobolev source func tion due to velocity field 
curvature was also fo und by Sellmaier et al.l ([1993), and 
lOwocki fc Pulsl ()1999[ ). We note that in the case of the res- 
onance lines plotted in Fig. [2] the line source function at 
large r radii is roughly propo rtional to S* ~ rij/ni ^ r~^ 
(e.g., iKudritzkifc Puls 200m. 

The minimum of c'~^^^ close to the star is also connected 
with the velocity gradient changing significantly within the 
resonance zone. Thus, a given line also picks up the radia- 
tion corresponding to a lower velocity gradient leading to a 
further reduction in the radiative force. For velocities com- 
parable to or higher than the ion thermal speed, the lines 
are deshadowed because of the Doppler effect, and the ra- 
diative force increases. We note that we only include the 
thermal broadening, hence these effects occur for velocities 
lower than the sound speed. 

As the wind accelerates, the ratio of the CMF to Sobolev 
line force increases and reaches a value close to one for ve- 
locities higher than the thermal speed of the wind driving 
ions (Fig. [T]). This is unsurprising, because the Sobolev ap- 
proximation is applicable to regions with a large velocity 
gradient, which exist already close to the sonic point Vr = a. 
Owing to line overlaps, c^^^ is less than one in the outer 
wind regions, where it reaches only 0.7 — 0.8. 

To test the influence of line overlaps, we calculated the 
radiative force with only 50 carefully selected optically thick 
lines that do not overlap (see Fig. [3]). The pronounced 
minimum for velocities lower than the sound speed is still 
present here, but in the outer regions the value of c'~^^^ is 
approximately one, supporting the validity of the Sobolev 
approximation for supersonic velocities. 

To understand more clearly the infiuence of line overlaps 
on the radiative force, we constructed another artificial line 
list using our set of non-overlapping lines. Each line in this 
new line list is counted twice with all parameters being 
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CO 




100 



Fig. 2. The line source function for resonance lines of O iv 
at 790 A, Sv at 786 A, and Fev at 388 A as a function of 
relative wind velocity. The source function is plotted rela- 
tive to its value at the sonic point Sa for the wind model 
of HD 93204. 



1.2 



0.6 



0.4 



with overlaps 
no line overlaps 



4 
r/R* - 1 



Fig. 3. The radial variation in the ratio of the CMF to 
Sobolev line forces in the wind model of Tc 1 central star 
with and without line overlaps. 



1.4 
1.3 
1.2 

1.1 
1 

0.9 
0.8 
0.7 
0.6 
0.5 



a / 



0.001 



0.01 0.1 



10 



Fig. 4. The ratio of the CMF to Sobolev line forces at ra- 
dius r = 5 i?* in the Tel wind model in the dependence 
on the line shift. The value of hydrogen thermal speed is 
denoted in the figure. 



completely the same, however with a line center shifted by 
Vijl^Vi/c, where Auf is a free parameter. For Aw^ ^ a, 
all twin lines completely overlap leading to a significant 
decrease in the radiative force with respect to the Sobolev 
one that does not account for the line overlaps (see Fig. 14]) . 
For Au^ > a, the lines at a given point do not overlap, 
but one of the twin lines "sees" the flux absorbed by the 



Table 2. Comparison of calculated wind parameters de- 
rived using CMF and Sobolev line forces 



Star 



HSF" 



M 
[Me year" 



[kms 



li 



CUV 



M 
[Mq year" 



[kms-^l 



First stars 



M500-1 


6.3 X 10""** 


2750 


6.1 X 10-"« 


4010 


M500-2 


4.0 X 10"°^ 


1930 


2.0 X lO"'^'^ 


1310 


M500-3 


2.1 X 10-°^ 


580 


1.4 X 10-°'^ 


790 


M500-4 


3.8 X 10"°® 


600 


2.8 X 10"''* 


620 


O stars 


HD 24912 


4.4 X 10"' 


2270 


2.3 X 10"' 


2030 


HD 37043 


6.2 X 10-^ 


2340 


4.1 X lO"'^ 


2000 


HD 47839 


2.2 X 10"^ 


3080 


1.0 X lO"'^ 


2970 


HD 54662 


7.9 X 10~^ 


2190 


4.1 X lO"'^ 


2050 


HD 93204 


1.3 X 10"^ 


2290 


5.9 X lO"'^ 


2080 


HD 149757 


4.7 X 10"* 


2040 


2.9 X 10"* 


1860 


HD 203064 


5.7 X 10-^ 


2080 


3.8 X 10"-^ 


1780 


HD 209975 


8.4 X 10"^ 


2430 


5.5 X lO"'^ 


1960 




Central stars of 


planetary nebulae 




NGC 2392 


3.7 X 10-« 


490 


1.8 X 10"* 


500 


NGC 3242 


3.1 X 10"'' 


2000 


2.0 X 10"-^ 


1890 


IC 4637 


3.1 X 10"* 


1440 


1.4 X 10~* 


1270 


IC 4593 


7.4 X 10"* 


730 


3.8 X 10"* 


660 


He 2-108 


9.5 X 10"* 


730 


4.7 X 10"* 


700 


IC418 


9.5 X 10"* 
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second line, leading to a reduction in the radiative force 
even in this case. For lS.vi « Woo, one of the lines is affected 
by the emission from the second one, leading to an increase 
in the CMF radiative force relative to the Sobolev one. 

A similar redu c tion i n the line force by multiline effects 
was found bv lPulsl ()1987D . We note that the multiline effects 
were also studied with respect to the multiple radiative mo- 
ment um deposition in Wolf-Rayet star winds (jGavlev et al.l 
|1995[ ). However, these effects are probably of minor impor- 
tance here due to the low density of the studied winds. 

The CMF radiative force, which is lower than the 
Sobolev one because of line overlaps, causes a decrease in 
the mass-loss rate of CMF models with respect to Sobolev 
ones (see Table [5]) . The ratio of CMF to Sobolev mass- loss 
rates is about 0.58. The only exception is the model M500- 
1, for which the CMF mass- loss rate is nearly the same as 
the Sobolev one. The reason is that the star is so hot, that 
the wind is accelerated mainly by a dozen Ov and Ovi 
lines. For a critical point velocity, these lines do not over- 
lap, hence c^^^ w 1, and the CMF and Sobolev mass- loss 
rates are nearly the same. 

The resulting wind parameters of the central stars of 
planetary nebul ae can be compared with those derived from 
observations bv iPauldrach et aTl (J2004 see Fig. [5]). There 
is reasonable agreement between the wind parameters pre- 



Pauldrach et al. 



Pauldrach et al. 



dieted by ourselves and those derived by 

(I20QJ, see Fig. [5]). The mass-loss rates of 

(2004) are on average a factor of about 1.6 higher than those 

derived by ourselves. This is most likely partly because of 

the simplifications included in our code, e.g., the neglect 

of continuum opacity sources, and partly by the different 

abundances adopted. 
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Fig. 5. Comparison of our derived mass- loss rates (up- 
per panel) and terminal velocities (lower panel) of the 
centr al stars of planetary nebulae with those derived by 
iPauldrachet al- (, 2004). 
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Fig. 6. The mass-loss rate of HD 209975 in the models with 
additional turbulent line broadening relative to the models 
with zero turbulent velocity. 



6. Models with base turbulence 

The existence of a region close to the stellar surface where 
the CMF line force is low compared to the Sobolev one (see 
Fig. [T]) is partly caused by the source function gradients at 
the wind base and partly by the Sobolev approximation 
not being applicable to the subsonic regions. The CMF line 
force increases at the moment when the line starts to ab- 
sorb the radiation that has not been absorbed yet, i.e., the 



radiation from the line wing. Because up to now we have 
assumed pure thermal line broadening, the velocity width 
of low CMF line force is of the order of the metallic thermal 
speed (which is roughly 0.13a in the case of iron). The wind 
mass-loss rate in our models is determined in the region of 
supersonic wind, close to the critical point where the wind 
velo city approaches the speed of radiative-acoustic waves 
([Abbot t 1980 , Feldmeier et al. 2008). Thus, the region of 
low CMF line force close to the star does not significantly 
affect the wind mass-loss rate. 

However, if the line broadening were larger (due to sur- 
face turbulence), then the region of low CMF line force 
could spread out to large velocities comparable to the tur- 
bulent one. When the turbulent velocity is comparable to 
the critical point velocity, below which the wind mass-loss 
rate is set, this could cause a significant decrease in the 
wind mass-loss rate. To test this, we calculated wind mod- 
els with additional line broadening, which we attributed to 
the turbulent one. In this case, the line profile width is given 
not by Eq. (g]), but by 



Aiy^j = 



2kT 



turb 



(15) 



where Wtmb is the adopted turbulent velocity. 

The results of numerical models indicate that with in- 
creasing turbulent broadening the velocity width of low 
CMF line forc e increases leading to a lower mass-loss rate 
(see Fig. [HI cf. lLucvir2007D . Hence, in the presence of turbu- 
lence the wind parameters may not depend only on the ba- 
sic stellar parameters (effective temperature, radius, mass) 
but also on the line turbulent broadening. Moreover, this 
effect can possibly be one of the reasons why the mass- 
loss rates derived from observational analyses that t ake the 
clumping into account (Bouret et al. 2003, Martin s et al.l 
l2005l ) are systematically lower than the predicted ones. 

For velocities higher than a few times the turbulent one, 
the Sobolev approximation should be applicable. At these 
high velocities, one expects that the line force becomes close 
to the Sobolev one. Because now the same force (as in the 
model with zero turbulent broadening) accelerates the wind 
of lower de nsity, one exp ects the terminal velocity to in- 
crease (e.g.. lGavlevll200l[ ). becoming much higher than the 
observed one. However, our models do not predict a signif- 
icant increase in the terminal velocity Woo, which is in the 
range 1900-2200 km s'^ for the wind models of HD 209975 
with different turbulent broadening. This is caused by the 
stronger blocking of stellar radiation by increased line over- 
laps mainly in the region with v^ "^ a. We note that lines 
broadened by turbulent motions are able to block the flux 
more efficiently than lines broadened purely thermally. 

Observational studies consider the turbu lence already 
present in t he photospheres of O star s (e.g.. iBouret et all 
I2003L 120051 iMartins et all |2004 I2005D with turbulent ve- 
locities of about 2 — 25kms~^. Macroturbulent veloci- 
ties in B s upergiants may be even higher, about 30 — 
lOOkms-i (iHo warth et"aLl Il997l iMarkova fc P"^ I2008D . 
Convective layers and s urface pulsational mot i ons are also 



expec ted theoretically (jCantiello et al.l l2()09l. lAerts et al 
2009 ). Turbulence can spread in the wind (JFeldmeier et al. 



19971 ). leading to a decrease in the wind mass-loss rate, as 



shown here. We also note that many O stars exhibit turbu- 
lent velocities in the range 10 — 20 kms^^, where we expect 
a high sensitivity of the predicted mass-loss rate to the tur- 
bulent velocity (see Fig. [B]). 
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Fig. 7. The dependence of the radial velocity on radius for 
solutions with different boundary densities (mass-loss rates) 
close to the stellar surface for the model star 500-1. Each 
consecutive model (from down to up) differs by a factor of 
1.5 in the mass-loss rate (thin lines). The thick line denotes 
the unique solution that smoothly intercepts the critical 
point. The dashed line denotes solution with the Sobolev 
line force. 

The basic results presented here will, in the future, be 
tested in more detail using models that also account for 
the continuum opacity and CMF line source function in a 
separate study. 

7. The solution topology 

In Fig. [71 we plot solutions with different base densities 
(mass- loss rates). In generall, with increasing base density 
the wind velocity increases until the density reaches a maxi- 
mum value. There is no solution that is smooth out to large 
radii for the densities higher than the maximum one. The 
solution with maximum density is very similar to the criti- 
cal solution of Sobolev models (see Fig. [7]). Moreover, there 
are many solutions that smoothly pass through the sonic 
point V = a for different mass-loss rates. 

This indicates that the critical point of non- 
Sobolev models is close to the CAK critical point 
(jCastor. Abbott fc Klein 1 11975[ ) and that the sonic point 
is not a point where the wind mass-loss rate is determined. 
The reason is that even in the non-Sobolev models the ra- 
diative force is not given locally by wind density and veloc- 
ity, but depends on the wind properties in a close neighbor- 
hood of a studied point. This dependence on the non-local 
properties at its limit approaches the Sobolev approxima- 
tion for very thin resonance layers (for very large velocity 
gradients) . 

8. Conclusions 

We have presented hot star wind models in which the radia- 
tive force is calculated using the solution of the comoving 
frame (CMF) radiative transfer equation. The wind models 
were calculated for three different groups of stellar param- 
eters (corresponding to evolved first stars, O stars, and the 
central stars of planetary nebulae) to compare the CMF 
and Sobolev radiative forces for a broader range of stellar 
parameters. 

The comparison of the CMF radiative force with an ap- 
proximate one calculated by assuming the Sobolev approxi- 
mation showed that the Sobolev line force is slightly higher 



due to the neglect of line overlaps. Thus, the mass-loss rate 
of wind models that include the Sobolev line force is on 
average a factor of about 1.7 higher than a more realistic 
one calculated using CMF wind models. However, we note 
that the simple Sobolev approximation applied here for ref- 
erence can be improved to account f or line overlaps and 
continuum absorption (Olson 1982, Hummer fc Rybickil 
11985, Puis fc Hummer 1988, Pavlakis fc Kylafis 1996 |). We 
emphasize that modern hot star wind models include 
line overlaps l'e.g..lVink et al.l[200li iPauldrach et al.ll200lL 
iGrafener fc Hamannll2005[l . 

Without line overlaps in the case of purely thermal line 
broadening, the Sobolev approximation provides a reliable 
estimate of the radiative force even close to the wind sonic 
point. The CAK model therefore provides a good approx- 
imation for a solution obtained with non-Sobolev transfer. 
Below the sonic point, the CMF line force is significantly 
lower than the Sobolev one partly because of the strong gra- 
dients in the source function. This does not influence the 
mass-loss rate, as the wind mass-loss rate is set in the su- 
personic part of the wind below the critical point. However, 
when additional line broadening is present (e.g., the turbu- 
lent one) then the region of low CMF line force may extend 
outwards to the regions where the mass-loss rate is set. 
This results in a significant decrease in the wind mass-loss 
rate. We note that this is not a shortcoming of the Sobolev 
approximation because the Sobolev approximation is appli- 
cable to velocities higher than the turbulent velocity in this 
case. 

The influence of turbulent line broadening may cause 
a dependence of the wind mass-loss rate on the atmo- 
spheric turbulent motions. Because theoretical models are 
not yet able to predict the atmospheric turbulent motions 
in hot stars in detail, we are unable to provide reliable 
wind mass-loss rate predictions until the theory of the at- 
mospheric turbulence develops consid erably (however see 
ICantiello et"aI1l2009LlAerts et al.ll2009[ ). Nowadays, hot star 
evolution seems to be a deterministic one depending only 
on the initial stellar parameters, i.e., mass, metallicity, and 
rotational rate. However, as the properties of atmospheric 
turbulent motions seem to be a non-trivial function of stel- 
lar parameters, the evolution of hot stars may become less 
deterministic, becoming instead dependent on free param- 
eters describing the role of surface turbulence. 
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Appendix A: The solution of CMF radiative transfer equation 

To calculate the radiative force, the iMihalas et al.l ()1975[ ) method for the solution of the CMF radiative transfer equation 
is modified in such a way that th e v variable is sp ecified on the spatial grid, and u on the intermediate one. 

Following the notation of Mih alas et alj ()1975l) . the depth index d increases inward, ri = i?out > ^2 > • • • > tnd — R*: 
where Rout is the radius of the outer model boundary. The impact parameters p are labeled in order of increasing size by 
index j ■, pi = < P2 < ■ ■ ■ < Pnc < ■ • • < Pnc+nd, where NC is the number of rays intersecting the core. Along each ray 



with impact parameter pj, we define grid in z and optical depth r, Zji 



[Rlnt-P])''' > ^12 



2\l/2 



> 



Zj-,Ni, , where NIj = ND for j < NC, and NI^- = ND + NC + 1 - j for NC < j < ND + NC, r^i = < rj2 < • • • < rj,Ni, • 
The frequencies are labeled by index fc in order of decreasing values, v\> vi> ■ ■ ■ > i^nf- 

We assume that v is specified on the depth grid, and u is specified at intermediate grid points labeled by d ± ^ . 
Suppressing the ray index j in the following, wc define on each ray pj 



Xk.d+i/2 = 2 \x{vk,Zd+i) + x{vk,Zd)\ , 

^Tk4+l/2 ~ Xk,d+l/2 {Zd — Zd+l) , 

^Tk,d = 2 (^■^fc,d+l/2 + Art^d^l/2) • 



The difference form of the system of the equations in Eq. Q is 

U{vk, Zd+1/2) - u{'^k,Zd-l/2) 



^Tkd 



= v{vk,Zd) + 



lk,d 



v{Vk,Zd+l) -v{Uk,Zd) 
^'''fe,d+l/2 

where d = 2, . . . , NI — 1, and 



U{l^k, Zd+1/2) - S{Vk,Zd+l/2) 



A;^fc_i/2 ^ '^k-i-'^k, 

ad+1/2 
lkA+1/2 



^Vk-i/2 
lkA+1/2 



^Vk- 



1/2 



v{vk,Zd) - v(yk-i,Zd)\ , 

[M(l^fe,^d+l/2) --"(i^fc-l, Zd+1/2)] 



^d+l/2Xfc,d+l/2 



2 
Md+1/2 



/3<i+l/2Md+l/2J 



lk,d 



OLd 



rdXk,d 
Solving Eq. (jA.4bp for uli/k, Zd+1/2), we obtain 

v{iyk,Zd+i) -v{i'k,Zd) 



{l^^il + PdA) 



uii^k, Zd+1/2) 



where 



(1 + Sk-i/2,d+l/2)^TkM+l/2 1 + '^fc-l/2,d+l/2 

7A;,d+l/2 



'^fc-l/2.d+l/2 , . SiVk, Zd+1/2) 

U[L'k-l, Zd+1/2) + 



1 + Sk-1 



/2,d+l/2 



h- 



l/2,d+l/2 



^Vk-1/2 

Substituting Eq. (|A.8I) into Eq. (IA.4ap . we derive a linear system of equations for v{vk, Zd) 



^Tk., 



v{vk,Zd^ 



Arfc_d+i/2(l + 4-1/2, d+1/2) 



{l + 5k-l/2M)v{vk,Zd) 



v{vk,Zd) 



^TkA+l/2{'^ + '^fe-l/2,d+l/2) ATi^^.d-l/2(l + 4-1/2, d-1/2) 

v{vk,Zd-l) 



1 



A'''fc,d-1/2(1 + 4-1/2, d-l/2)_ 

S{yk, Zd+1/2) S{vk,Zd-i/2) 



^Tk.d V 1 + <5, 



3fc-l/2,d+l/2 



5fc-l/2,d-l/2 



- &k-l/2,dV{vk^l,Zd) + 



1 



A-Tfe,, 



4-1/2, d-l/2w('^fe-l, Zd-1/2) 4-1/2, d+l/2w(l'fe-l, 2:^+1/2) 



1 + 4-1/2, d-1/2 



1 + 4-1/2, d+1/2 



(A.1) 
(A.2) 
(A.3) 

(A.4a) 
(A.4b) 

(A.5) 
(A.6) 

(A.7) 

(A.8) 
(A.9) 



(A.IO) 



This system should be supplemented by equations corresponding to the boundary and initial conditions. At the outer 
spatial boundary Zout, we assume no infalling radiation, consequently u — v and we derive from Eq. (|7b[) 



x{v, r) dz 
or, in a difference form 

v{vk,Zd+i) -v{vk,Zd) 



1 dv{iy,p,z) dv{v,p,z) ^f . , ^ 

- lii^iP, z) T. = t>{v, r) - viv.p, z), 



dv 



^Tk. 



d+1/2 



= v{l'k,Zd){l + 4-l/2,d) - h-l/2,dVi'^k-l,Zd) - S{vk,Zd). 



(A.11) 



(A.12) 
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The infalhng radiation at the inner boundary is taken from the model atmosphe res. The initial solution for vi is derive d 
using the solution of the radiative transfer equation neglecting the velocity fields (jMihalas h Hummer| ll974'. Kubat 199%. 

The velocity derivatives at the grid points are approximated as in the hydrodynamical code (see Krticka & Kubia 

I2OOII Eq. (A. 4a) therein) . The derivatives in the middle points between grid points are calculated as the average of the 
derivatives at the grid points. 

The system of algebraic equations Eq. (|A.10p with bound ary conditions is solved using the LAPACK package 
(http: //www. cs . Colorado. edu/~lapack^ Anderson et al. 119991 ). 



